A gamma-ray signature of energetic sources of cosmic-ray nuclei 
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Abstract 

Astrophysical sources of nuclei are expected to produce a broad spectrum of isotopes, many of 
which are unstable. An unstable nucleus can beta-decay outside the source into a single-electron 
ion. Heavy one-electron ions, thus formed, can be excited in their interactions with cosmic mi- 
crowave background photons, in which case they relax to the ground state with the emission of 
a gamma ray. Repetitive cycles of excitation and gamma-ray emission can produce an observable 
feature in the gamma-ray spectrum with a maximum around 8 GeV (for iron). We find that the 
observed spectrum of Centaurus A is consistent with a substantial flux of nuclei accelerated to 
0.1 EeV. A characteristic 5-10 GeV (iron) shoulder in the gamma-ray spectra of various sources 
can help identify astrophysical accelerators of nuclei or set upper limits on nuclear acceleration. 

PACS numbers: 98.70.Sa, 95.85.Pw, 98.70.Rz, 95.85.Ry 
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Heavy nuclei are detected in high-energy cosmic rays although the evidence is still 
controversial at the highest end of the spectrum [2j. The sources are difficult to identify 
because the arrival directions are altered by the galactic magnetic fields. It would be very 
desirable to identify astrophysical nuclear accelerators by some other means, and some sig- 
natures in gamma rays and neutrinos have been proposed [3| . We describe a new signature 
that can be used to identify powerful sources of cosmic-ray nuclei: a characteristic "shoulder" 
in the spectrum of gamma rays. 

Astrophysical sources, such as jets of active galactic nuclei (AGN), produce cosmic rays 
in an environment of high photon density, which can cause some accelerated nuclei to dis- 
integrate. However, a non- negligible fraction of nuclei can escape [4j. Gamma-ray data 
indicate that acceleration of cosmic rays does take place in AGNs: line-of-sight interactions 
of cosmic rays provide a natural explanation for the observed spectra of distant blazars j^J, 
assuming the intergalactic magnetic fields are relatively small jg|. Models allow the photon 
density in the acceleration region to vary from a low enough for the nuclei to escape to a 
high enough for photodisintegration (4]. It is reasonable to assume that nuclei coming out 
of the source are subject to some degree of photodisintegration, which would render them 
generally unstable. On average, it takes one beta decay for an unstable nucleus to become 
stable. 

Astrophysical sources can produce a number of different isotopes. Unstable nuclei decay, 
mostly, via beta decays. We will use iron as an example. For a fully ionized large nucleus, 
the beta-decay electron is created in a bound state with probability that is of the order of 
one 7] • Hence, an order-one fraction of nuclei turn into one-electron ions at distance CTd eC ay7 
from the source, where 7 = E/M, E is the nucleus energy, and M is its mass. Once a beta 
electron is bound in an ion, it is likely to remain bound even in the event of a subsequent 
beta decay, because the nuclear recoil and the relative change in the charge Z of the nucleus 
are both small. 

The relevant energy corresponds to 7 ~ 10 6 — 10 7 , for which the cosmic microwave 
background (CMB) photons excite the electron from its ground state to one of the higher 
bound states. The CMB photon energy in the frame of reference of the ion is ~ 7 keV 
for iron. We will see that the ion is ionized, on average, about one year or more after it 
was created in a bound beta decay. De-excitation of the excited state, which is practically 
instantaneous, produces a photon with energy E 1 ^ m 7 keV in the frame of the ion, which 
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corresponds to energy i? 7) iab ~ in the laboratory (CMB) frame. Repeated cycles of 

excitation and emission can create a gamma-ray signal that can help identify the sources 
of UHECR nuclei. We note that the probability of a beta electron to be captured on an 
excited nS atomic level (with n > 2) is approximately 0.2 of the capture probability for the 
ground state (n = 1) [7]. The de-excitation of these higher levels thus creates approximately 
0.16 gamma photons per ion at the time of beta decay, before the re-excitation of the ion in 
collisions with the CMB. 

The observable flux of gamma-rays emerging from the described process of excitation 
and de-excitation of the one-electron ions depends on the excitation rate T e of the ion from 
the ground state to a discrete excited level due to interactions with CMB photons, and the 
rate of the photo-ionization due to interactions with CMB photons (T;) or extragalactic 
background light (EBL) photons (T^ ebl) • The former process generates the gamma-ray 
flux, while the rate of ionization determines the lifetime of the ions in the photon bath. In 
what follows we refer to the laboratory (CMB) frame values of the rates. The rates can 
be calculated using the dipole approximation for the radiative transitions, since at Z w 26 
(iron) the atomic dynamics is still sufficiently nonrelativistic. In this approximation only the 
excitation of the ground state to P-wave states, discrete and in the continuum, is relevant. 
For the transition from the ground IS* to a specific nP state the rate in a photon bath can 
be readily found in the CMB frame in the form 
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In this expression u and k are the energy and the momentum of the photon in the bath and 9 
is the angle between the momenta of the photon and the ion in the CMB frame as well as the 
photon thermal distribution function 7Zcmb(^) = [exp(w/T) — at T = Tcmb = 2.725 K, 
while the excitation energy A n = e(nP) — e{lS) and the dipole matrix element (nP\f\lS) 
refer to the rest frame of the ion, and the sum runs over the polarizations of the nP state. 
Finally 7 = E/M is the 7 factor for the ion in the CMB frame. After angular integration 
the expression (JTJ) takes the simple form: 
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where 

X (z) = -z 2 In (l - e~ z ) . (3) 

It can be noticed that the lower limit in the u integral is the lowest photon energy at which 
the excitation to the energy A n is kinematically possible (in a head-on collision). 

For a simple hydrogen-like one-electron ion the dipole matrix elements in Eq.([2]) for 
discrete bound states can be readily found in terms of the rescaled Bohr radius az = clb/Z = 
(rrieaZ)' 1 , using the explicit form of the hydrogen wave functions: 

El(nP|r1lS)| J = ^ (^) 2n 4, (n>2) (4) 

with the excitation energy, clearly, expressed through the re-scaled Rydberg energy ez = 
Z 2 Ry = m 2 e a 2 Z 2 /2 as A n = (1 - n~ 2 ) e z . 

Using Eq.(|2]), one can find the total excitation rate T c = 1/A C of the ion to bound states: 

r e = ^a%T>x (2.) = 1 7 x (2.) ni\ 2 , (5 ) 
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where 

X W = 256 g^(^l) 2 ",[(l-0/ r ] , (6) 

with 7^ = ez/(2T) = 1.96 x 10 (Z/26) 2 being the characteristic value of the gamma factor in 
the present discussion. Function X(r), computed numerically is shown in Fig. 1. The series 
in Eq. ([6]) converges quickly, and even the first term provides an excellent approximation to 
the sum. 

Excitation of the states higher than 2P can result in the atomic cascade emission of a 
gamma ray with energy boosted from at least A 2 = (3/4) ez- However the contribution of 
the higher excitations to the rate T e is small. We, therefore, use the approximation that 
every excitation results in emission of a photon with energy (3/4) ez = 6.90 (Z/26) 2 keV, 
which, boosted to the laboratory frame, produces a uniform spectrum of gamma rays with 
energies up to (3/4) 7 e z = 101 GeV (Z/26) 4 (7/7^)- 

The rate of the ionization process in which the 'active' ions are destroyed can be calculated 
similarly to the discrete excitation, in terms of the atomic ionization function 
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FIG. 1: The functions X(r) and Y(r) describing the rates of excitation and ionization, respectively, 
versus the ratio r = 7/7^. 

where \p) is the state in continuum in the Coulomb potential whose momentum at infinity 
is asymptotically p. This function is well known in the theory of hydrogen atom, and can be 
read off from , e.g. Ref. [8j: S(A) = a 2 z e z s(A/e z ), where the dimensionless function s(w) 
is given by 
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and where the branch of the arctangent must be chosen to take the values between and n. 
Similar to Eq.(j2J), the ionization rate F{ in the CMB thermal bath, including the ionization 
to all energies above the threshold can be written as 
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y(r)=r s W ,P^. (10) 
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The ionization function Y(r) is readily calculated numerically, and the plot is shown in 
Fig. 1. Clearly, the ratio of the excitation and ionization factors determines the average 
number of gamma photons emitted by an ion during its life time Tj = 77 = X/Y. The 
plot for this ratio is shown in Fig. 2. 



FIG. 2: Two of the functions that determine the average number of gamma photons generated 
during the lifetime of an ion, iV 7 = minjry, fj, fj}. 



In addition to interactions with CMB, ionization can occur due to interactions with extra- 
galactic background light (EBL), with energies between 10 _3 eV and 1 eV in the laboratory 
frame. Since all of these photons have relatively high energies, they ionize the ions and termi- 
nate the cycle of absorptions and emissions. The cross section for the relevant energy range 
is a ~ <tt = 0.6 barn, and the number density can be integrated numerically using Ref. |9|. 



There are non-negligible uncertainties in the level of 



in the range (0.5 — 1.1) cm~ 3 , depending on a model For a mid-range reference value of 
0.8 cm -3 , the ionization by EBL photons occurs at a rate r^. ebl ~ V\ebl ~ 1/(0.7 Mpc). 
For a nucleus with a gamma factor rjz that passes distance d (which can be the distance 
to the source, or the length of a diffusive path in turbulent magnetic field), the number of 
photons emitted in its lifetime cannot exceed r\ = X/Y, or fj = Xi-ebl/K, or fj = d/X e , so 
the total number of photons per nucleus is 

N y = mm{r],fj,fj}. (11) 

As shown in Fig. 2, for a sufficiently distant source, the number of photons is maximized 
around r = 0.055, which corresponds to the photon energy around 8 GeV. 

According to Eq.flS}, the ion lifetime 7* depends on its gamma factor. Using the calculated 
function Y (Fig. 1), we estimate that this lifetime is minimal at 7 ~ 0.7757^, where Y = 
F max w 0.457, which for Z = 26 corresponds to r, ~ 5.4 x 10 7 s ~ 1.7yr and 7 ~ 1.14 x 
10 7 , which gives the upper bound of the gamma photon spectrum at E 7 ~ 79 GeV. The 
phenomenological significance of the discussed effect depends on the relation between the 
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lifetime of the ion 73 and the duration tn of the process in which heavy nuclei are generated 
and destroyed or decelerated. If tat ^> Tj, the flux of the ions Fi and the flux of the gamma 
photons F 7 are related as 

F 1 = N n F i (12) 

In the opposite case, *< t,, the rate of emission of the gamma photons exponentially 
decays and their flux is given by F 1 = T e e~ Fit J Fi dt, where the integral is over the duration 
of the burst flux of ions. However, for the mechanism in question, the time scale for creation 
of the one-electron ions is determined by the beta decay times rp ~ 1 4- 10 4 s, dilated by the 
factor 7, which results in the effective emission times for the ions that are longer than Tj. 
Hence, in our subsequent estimates of the gamma-ray flux we will use Eq. (|12p . 

We now compare our predictions with the data for diffuse background and for point 
sources. Diffuse background is probably dominated by the nuclei diffusing in the magnetic 
field of our Milky Way galaxy. The flux of nuclei at the relevant energy ~ 10 17 eV can 



be inferred from the data 
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F Fe = 2 x lcr^cnrVV -1 . (13) 



This implies the local density of nuclei np e = 1-3 x 10 -23 cm -3 . To translate this into 
the photon flux, one has to make assumptions about the largely unknown geometry of 
the diffusive magnetic fields in our galaxy. The simplest estimate, assuming a spherical 
distribution of diffusing nuclei of the size L, in which each nucleus spends time r yields 
diffuse photon flux at energies (5-10) GeV 
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This is well below the diffuse background reported by Fermi 

We now discuss the predictions for point sources. One particularly interesting candidate 
is Centaurus A (NGC 5128), which is the closest active galaxy at distance d = 3.4 Mpc. A 



cluster of UHECR events in PAO data appears to favor Cen A as the origin of UHECR 11 ]. 



At lower energies, the galactic magnetic fields isotropize the arrival directions. It is pos- 



sible that all or most of UHECR come from a single nearby source 
sources 



12j, or from galactic 



131 ] . As an example, we apply our discussion to Cen A. The observed flux in Eq.(iT3~ 



can be used as an upper limit on the Cen A contribution. If the 5-10 GeV shoulder that 
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FIG. 3: Signature of nuclear emission for Cen A (solid line) and the data from Fermi 
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HESS [15I . [la ]. The flux normalization corresponds to the observed flux Fp e (E > 10 Ll cV) — 
2 x 10 _14 s _1 cm _2 sr _1 Based on this comparison alone, the upper limit on Cen A luminosity 
in nuclei with energies above 10 17 eV is 4 x 10 42 erg/s, close to the observational upper limit. 



we have described was ruled out by the data, one could set a stronger upper limit on nu- 
clear acceleration in Cen A. The emission of gamma rays is possible in a narrow band of 
energies, and the flux is enhanced around r = 0.055, as shown in Fig. 2. The dominant 
gamma-ray energy is r7^(7keV) = 8 GeV. In this case the expected flux of 8 GeV photons 
is .FcenA ~ 47T Fp e Nj, where we take A^ 7 = 10 3 , consistent with Fig. 2. This flux is in agree- 
ment with the data, as shown in Fig. 3. The spectral feature at 5-10 GeV may, in fact, 
result from the effect we have described. However, there are good reasons to believe that 
local, Milky Way sources are responsible for most of the nuclear flux in this energy range 

Using the data from Fermi and other gamma-ray instruments, one may be able to identify 
powerful nuclear sources at large distances. The shape of the "shoulder" due to cosmolumi- 
nescence may contain information about the range of isotopes among the nuclei. 

In summary, we have described a new signature of astrophysical nuclear accelerators, 
which can be used to identify or constrain astrophysical sources of nuclei. 
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